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^ ; Abstract 
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Based on a counting scheme that exphcitly takes into account the large momentum 
^ , ^/Mrn^ characteristic for pion production in nucleon-nucleon collisions we calculate all 

' diagrams for the reaction A^A^ — > NNtt at threshold up to next-to-leading order. At this 

o 

00 
O 
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' order there are no free parameters and the size of the next-to-leading order contributions 

is in line with the expectation from power counting. The sum of loop corrections at that 
I order vanishes for the process pp ppn^ at threshold. The total contribution at next- 

' to-leading order from loop diagrams that include the delta degree of freedom vanishes at 

threshold in both reaction channels pp — > ppn^, pn7r~^. 
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I The high precision data for the processes pp ppir^, pp pmr'^ and pp — dir^ in the 

threshold region [|l| have spurred a flurry of theoretical investigations. The first data on neutral 
pion production were a big surprise because the experimental cross sections turned out to be 
a factor of five larger than the theoretical predictions based on direct pion production and 
neutral pion rescattering fixed from on-shell nN data 0, Subsequently, it was argued that 
heavy-meson exchanges might be able to remove this discrepancy 0. On the other hand, it 
was found ^ that the (model-dependent) off-shell behavior of the full tcN T-matrix can also 
enhance the cross sections near threshold considerably. 

Due to their nature as pseudo-Goldstone bosons the dynamics of pions is largely constrained 
by chiral symmetry. Thus one might hope that effective field theory studies which incorporate 
these constraints strictly will help to resolve the so far confusing situation. In the literature there 
are several calculations carried out in the framework of tree-level chiral perturbation theory 
including the dimension two (single-nucleon) operators for neutral pion production |0, |], |, [10| 
as well as for charged pion production [|IT], . A common feature of these calculations is that 
the contributions from the isoscalar pion rescattering interfere destructively with the direct 
production amplitude, thus leading to an even more severe discrepancy between experiment 
and theory. It should be noted that such an interference pattern is in contradiction to the one 
found in phenomenological approaches P, H. Furthermore, within the Weinberg scheme, where 
all momenta are considered of the order of m^, one loop calculations have been performed 
for neutral pion production pp — >■ ppir^ [|^, According to some of these works the 



loop corrections are larger by at least a factor of two compared to the tree level diagrams, that 
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according to the counting scheme apphed appear one order down. This feature (if correct) would 
seriously question the convergence of the chiral expansion for pion production in NN-coUisions. 



On the other hand, according to ref . |T6[ the chiral expansion seems to show convergence in the 
case of p-wave pion production. 

The purpose of the present work is to present a complete next-to-leading order calculation 
of the reaction A^A^ — > NNn at threshold. In particular, we evaluate all one-loop diagrams at 
next-to-leading order employing a counting scheme that takes into account the large momentum 
yJMrriT, characteristic for pion production in NN-collisions, as suggested in refs. |T^. We 
consider also the contributions from explicit delta-isobars at tree level and at one-loop order. 
To the order we are working there are no free parameters and we demonstrate that the size of 
the individual next-to-leading order contributions is in line with the expectations from power 
counting. 

Let us begin with writing down the general form of the threshold T-matrix for the pion 
production reaction Ni{p) + N2{—p) —>■ N + N + n in the center-of-mass frame, which reads 



m 



B 



with (Ti 2 and ri 2 the spin and isospin operators of the two nucleons. denotes the three- 
component isospin wave function of the final state pion produced in an s-wave state, e.g. 
= (0, 0, 1) for 7r°-production and (p = (1, i, 0) / a/2 for vr """-production. The complex amplitudes 
A and B belong to the transitions ^Pq ^So and ^Pi ^Si in the two-nucleon system, 
respectively. In fact the selection rules which follow from the conservation of parity, angular 
momentum and isospin allow only for these two transitions for the reaction A^A^ NNtt at 
threshold. In the case of neutral pion production pp ppn^ the threshold amplitude A is the 
only relevant one whereas in charged pion production pp — >■ pmr^ both threshold amplitudes 
A and B can contribute. Note that the threshold T-matrix written in eq.(l) incorporates the 
Pauli exclusion principle since combined left multiplication with the spin-exchange operator 
(1 + ai ■ (T2)/2 and the isospin exchange operator (1 + ri ■ r2)/2 reproduces T^f^{NN — > NNn) 
up to an important minus sign. The magnitude of the nucleon center-of-mass momentum p 
necessary to produce a pion at rest is given by: 



\P\ 



^m^iM + mjA), (2) 



with M = 939 MeV and m,r = 139.6 MeV denoting the nucleon and pion mass, respectively. 
Eq.(2) exhibits the important feature of the reaction A^A^ — > NNn, namely the large momentum 
mismatch between the initial and the final nucleon-nucleon state. This leads to a large invariant 
(squared) momentum transfer t = —Mm.,^ between in- and outgoing nucleons. The appearance 
of the large momentum scale y/Mrn^ in pion production demands for a change in the chiral 
power counting rules, as pointed out already in ref.0. In addition, it seems compulsory to 
include the delta-isobar as an explicit degree of freedom, since the delta-nucleon mass difference 
A = 293 MeV is comparable to the external momentum p ^ \fMrn^ = 362 MeV. The hierarchy 
of scales 

M > p ~ A > , (3) 

suggested by this feature is in line with findings within meson exchange models where the 
delta-isobar gives significant contributions even close to the threshold [P 
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Let us now state our counting rules. The external momentum p ~ y/Mm^ sets the overall 
scale relevant for the process A^A^ — > NNtt. This momentum scale p enters the internal lines 
of tree and loop diagrams. Therefore we count all four- moment aP^ inside loops generically 



as order p and the loop integration measure / dH as order j9^. A pion propagator is counted 
as order The delta-propagator of the form 1/ {energy — A) counts as order 1/p, since we 

made the choice A ~ p. For the nucleon propagator of the form 1/ energy one has to distinguish 
whether it occurs outside or inside a loop. The associated residual energy counts as order 
outside a loop and as order p ~ y/Mm^ inside a loop. Furthermore, external pion energies are 
counted as order m^. 

According to these counting rules one-loop diagrams contribute at order p"^ in the expansion 
of the T-matrix and thus generate threshold amplitudes of the form S ~ p ~ y/M^. The 
new counting rules demand also for a reordering of the terms in the interaction Lagrangian, 
since " relativistic corrections" proportional to nucleon kinetic energies p'^/M are now of the 
same order as "leading order contributions" proportional to residual nucleon energies. Several 
examples of this effect will be encountered here. 

In Fig. 1, we display tree-level diagrams which according to the abovementioned counting 
rules contribute at leading order, next-to-leading order and next-to-next-to-leading order. Dia- 
grams for which the role of both nucleons is interchanged and diagrams with crossed outgoing 
nucleon lines are not shown. Subsets of four diagrams obtained by these operations map 
properly onto the crossing antisymmetric threshold T-matrix eq.(l). Diagram a) involving the 
(isovector) Weinberg- Tomozawa vrTrA^A^-contact vertex gives a leading-order contribution of the 
form: 

with gA — 1.3 the nucleon axial vector coupling and /tt = 92.4 MeV the pion decay constant. 
It is important to note that the Weinberg- Tomozawa vertex generates here a proportionality 
factor at "leading order" in the chiral vrA^-Lagrangian via the pion and nucleon (residual) 
energies as well as through a "relativistic correction" of the form p^ /M. This factor of m,r gets 
finally canceled by the pion propagator [mT^{M + m,r)]~^- Obviously, the isovector Weinberg- 
Tomozawa vertex cannot contribute to the neutral pion production threshold amplitude A. 
From the one-pion exchange diagram b) one finds: 

This result stems from the recoil correction to the irNN-veitex proportional to {rriT^/M) Bi ■ p 
with the m^-factor getting now canceled by the intermediate nucleon propagator. Furthermore, 
the product of the two vertices on the left nucleon line (ai -pY = Mm^^ is canceled by the pion 
propagator. The ratio B'^^'^^ / A^^^'^ = 3 has its origin in the isospin factor of diagram b). From 
the analogous diagram d) with one virtual delta-isobar excitation one finds: 

. (A) _ 9Am^ g(A) _ Q (Q. 

4M/3A' « U, W 

where we have used the empirically well satisfied relation Ha = 3gA/V^ for the vrA^A-coupling 
constant. The spin and isospin transition operators entering the vrA^A-vertex (/i^/2/^)S' ■ pTa 

"i^^Baryon energies are residual energies with the nucleon mass M subtracted. 
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satisfy the usual relations SiSj = {26ij — itijkCTk)/^ and T^tJ = (26ab — i^abc'^c)/'^- The latter 
isospin relation is the reason behind the vanishing of B^^\ According to our counting of the 
mass-splitting A the term A^^^ in eq.(6) is a next-to- leading order contribution, since A ~ p (c.f. 
relation (j^)). Diagram f) involves the second order chiral TrvrA^A^-contact vertex proportional 



to the low-energy constants Ci^2,3,4 ^^e following contributions to the threshold 



amplitudes at next-to-next-to-leading order: 

-4'=-' = ||7|fe + 2c,-4c0. e<'.. = ||^(c. + C3 + 2c,-4cO. (7) 

In a previous calculation in ref.0 (see eq.(32) therein) the C2-term has been found with a relative 
factor 1/2 smaller. The reason for this discrepancy is again that "relativistic corrections" from 
the C2- vertex are of the same order as its "static" contribution, since p'^/M = m^r. We also note 
that our results eqs.(4-7) agree up to the respective order with those of the fully relativistic 
calculation in ref.|T^ where no approximations to the threshold kinematics have been made. 



We do not specify here the contributions from diagrams c), e) and g) in Fig. 1 which are 
proportional to the (a priori unknown) strengths of four-nucleon contact-vertices etc. It is 
important to note that already at leading order long-range effects from pion-exchange and 
short-range contributions appear simultaneously. 

Let us now turn to the non-vanishing one-loop diagrams at threshold. Not every loop 
diagram appearing formally at next-to-leading order truly contributes at that order. In case of 
the diagrams b) and c) in Fig. 2 the (spin-independent) one-loop vrA^-scattering subdiagrams 
are proportional to m^, and this pushes their contributions to the threshold T-matrix eq.(l) 
beyond next-to-leading order. A closer inspection of diagrams a) in Fig. 2 reveals that they 
contribute in the form m-^lnrnT^ to the threshold amplitudes A and B, i.e. beyond next-to- 
leading order. The specific vertex structures of diagrams d) and e) in Fig. 2 make also their 
next-to-leading order contributions vanishing. Therefore we have to focus only on the diagrams 
shown in Figs. 3 and 4. 

We evaluate only the genuine next-to-leading order pieces of the loop integrals emerging 
from the diagrams in Figs. 3 and 4. For instance, in the integrands we can systematically drop 
terms of order compared to Iq (and A). Straightforward but tedious evaluation leads to the 
following next-to-leading order contributions of the one-loop diagrams in Fig. 3 with nucleons 
only: 



^iN-ioo,) ^ i^^^{-2 - 1 + 3) , = ^^^P^(-2 + + 3) . (8) 

Here, the numerical entries correspond to the diagrams a), b) and c), in that order. Interestingly, 
the total next-to-leading order loop contribution vanishes identically for neutral pion production 
grams a) and c) in Fig. 3 have been calculated fully relativistically (i.e. 
without any approximation to the threshold kinematics) for pp — > ppir^ in ref.[0. It is an 
important check for our calculation that the non-analytical piece proportional to ^/Wm^ agrees 
with the one derived by expanding eq.(16) in ref. |T^. In addition, after correcting a sign error 



in ref . |jTj] and extracting at threshold the truly next-to- leading order pieces from that work our 



results agree with theirs [pi 



Numerically, the loop correction in eq.(|D gives 

j^{N-ioop) ^ g\./M^/2hQfl ~ 0.70 fm^ 
This is about 50% of the leading order one-pion exchange contributions ~ 1.33 fm^ or 

^(it) ~ 1.69 fm*^. Indeed from chiral power counting one expects a similar suppression factor 
p/M = ^m^/M ~ 0.4. 
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According to our counting of the mass difference A ~ \fMrn^ loop diagrams with exphcit 
delta-isobars are of the same order as those with nucleons only, namely of order p^. The relevant 
one-loop diagrams which generate truly next-to-leading order contributions are shown in Fig. 4. 
Straightforward but tedious evaluation leads to the following result: 

^(A-Zoop) ^ ^|||^(8 - 12 + 1 + 3) , ^(A-^oop) ^ ^|||^(8 -12 + 3 + 1) (9) 

with the numerical entries corresponding to the subclasses a), b), c) and d), in that order. The 
relevant combination of loop functions reads: 

ir(A) = 2 Jo(-A) - 2A /o(-Mm,) + (2^^ - Mm^) 7o(-A, -Mm,) , (10) 



with the following loop integrals |20| truncated at lowest order according to our counting scheme: 



Jo(-A)=4AL(A) + A(^ln^-i) ^O(p), (11) 

1 / Mm \ 

loi-Mm^) = -2L(A) - ^ (^1 + In ^ j ~ 0(/) , (12) 



70 -A, -Mm, = , \ ——^ arctan % ~ 0{j>~^) . 13 

^TX^^Mm.,, Jo 1 + 2A 



The (scale dependent) quantity: 



A 



d-4 



167r2Lrf-4 2 



+ -(7b - 1 -ln47r) 



(14) 



denotes for the standard divergent piece in dimensional regularization. The reason for grouping 
together the three specific diagrams into subclass b) is that this way the (in the chiral limit) 
singular term A^Jo(— A)/Mm, does not appear explicitly. Evidently, the formal limit A — 
corresponds to loops diagrams with nucleons only, and therefore K{0) = — v^Mm^/16 enters 
eq.(8). Note, however, that for planar box diagrams this limit becomes inconsistent with the 
counting scheme employed. 

One concludes that the contributions stemming from loop diagrams with delta-excitation 
vanish identically for both threshold amplitudes A and B, respectively for both reaction chan- 
nels pp —>■ ppiT^ , pnTc^ . The complete cancellations in eq.(9) are actually important consistency 
checks for our power counting scheme A ~ p. The combination of loop functions K{A) ~ p 
in eq.(lO) is divergent, but at next-to-leading order there is no local counter term to absorb 
divergences. 

In summary, we have performed here a complete next-to-leading order calculation of the 
reaction A^A^ NNtt at threshold. We have employed the counting scheme developed in refs. 
[§, |TB|, that explicitly accounts for the large momentum p ~ y/Mm^ characteristic for this 



process. We find that the total next-to-leading order loop corrections either vanish or are in 
accordance with the expectation from power counting. At this stage we conclude, that the chiral 
expansion seems to converge also in the s-wave. Note, however, that at next-to-next-to-leading 
order a large number of loops enters, that have not yet been evaluated completely. 

In order to compare our results directly to pion production data the emerging chiral operators 
have to be folded with (realistic) NN-wave functions. This convolution has been carried out 
in ref.|]T3[ in a way that the symmetries are preserved. However, in that work the traditional 
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Weinberg counting has been used. Consequently the results presented in ref. ||T5[ do not allow 
any firm conclusion about the convergence of the chiral series, since contributions of different 
orders are mixed and the next-to-next-to-leading order is incomplete. Based on our counting 
scheme a complete next-to-next-to-leading order calculation is within reach and should be 
performed. Another direction should be the calculation of loop corrections to the higher partial 
waves amplitudes. 
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Figure 1: Tree level contributions to threshold pion production at leading order (a,b,c), next- 
to-leading order (d,e) and next-to-next-to-leading order (f,g). A single solid, double solid and 
dashed line denotes a nucleon, delta-isobar and pion, respectively. Leading (subleading) order 
vertices are symbolized by solid dots (open circles). 
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a) 



c) 



b) 



d) 



e) 



f) 

Figure 2: One- loop diagrams that start to contribute at next-to-next-to-leading order. For 
further notations see Fig. 1. 




Figure 3: Ncxt-to- leading order one-loop diagrams for pion production at threshold with nu- 
cleons only. For further notations see Fig. 1. 
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a) 



b) 




Figure 4: Ncxt-to-lcading order one-loop diagrams for pion production at threshold with inter- 
mediate delta-isobars. For further notations see Fig. 1. 
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